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• RELATION OF "LILIENTHAL EFFECT" TO DYNAMIC SOARING FLIGHT.* 

By Rodericti Fick. 

Otto Lilientlial observed that a flat surface, capable of ro- 
tation around a horizontal axis and supported by a counterpoise, 
oscillates vertically in a natural wind, whereby the momentary 
mean value of all the positions of said surface always forms an 
angle of 3-4° above the horizontal. The attempted explanations 
of this "upward component" long remained fruitless and finally 
led to doubting the results themselves and regarding them as er- 
rors of observation. On the assumption, however,, that the ob- 
servations of the upward component were correct, the turbulence 
of the wind was the only possible cause of the phenomenon. Nei- 
ther the Knoller-Betz theory nor the variations of the wind in 
horizontal flow explained the effect, so long as flat or symmet- 
rical profiles were employed, but, even in these cases, the ef- 
fect was said to be always produced. 

In this article, the phenomenon of the upward component in 
the case of flat surfaces will be referred. to as the "Lilienthal 
effect .J* Otto Lilienthal' s further observation that a cambered 
surface,, fixed at a suitable angle, experiences a forward thrust 
toward the wind, I will l^ere distinguish : from the "Lilienthal 

e ffect" by calling it the " Hnoller-Bet z effect," because Knoller 

* From "Zeitschrift fur Flugtechnik und Motorluf tschiff ahrt , " 
1924, November 28, pp. 244-246, . and December 12, pp. 258-260. 
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and Bet 2 first satisfactorily explained this effect hy means of 
vertical wind oscillations, at the basis of which, however, there 
lie air motions which the Lilienthal effect could not produce- 
Th© Lilienthal effect is also produced in the case of cambered 
airfoils. It will be shown that the same cause, which can pro- 
duce the Lilienthal effect on flat airfoils, considerably 
strengthens the thrust due to the Knoller-Betz effect- 

R. Knoller and Schmauck have already discovered the condi- 
tions of motion of the air, which can explain the Lilienthal 
effect (Flug- und Motortechnik, Vol. Ill, No. 22, "Die Gesetze 
des Luftwider standee," by R. Knoller; "Beitrag zur Erklarung des 
Segelf luges" by Schmauck, "Z5eitschrif t fur Flugtechnik und Mot.or- 
luftschiffahrt 1913). ETeither one, however, discovered absolute 
continuity., but only approximate or limited. Knoller assumes un- 
symmetrical waves, in which the downward phases are steeper tlian 
the upward phases, while the latter are correspondingly longer, 
the flow cross-section and velocity being constant. Schmauck, in 
an extension of the work of Betz, takes into consideration, in 
addition to the oscillations in direction, the simultaneous os- 
cillations in strength- and thereby emphasizes the special case, 
in which the oscillations in direction and strength coincide in 
such manner that the upward flow is faster and the downward flow 
slower- The question now arises as to whether there are meteoro- 
logical causes for such motion conditions of the air. 

Schnell, of Munich, has proposed a method for calculating 
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ttie energy of tur/bulsnce, necessarily existing near tlie ground, 
f rom tlie increase in the wind velocity with increase in height 
above the ground, "by friction on ground obstacles (Lecture "before 
the Aviation Convention' at Munich in the fall of 1923). Accord- 
ing to Schnell, the direction of transmission of the energy of 
turbulence in the higher layers of air is necessarily vertical 
and consequently exerts a lift on every object in the air. The 
transmission itself does not show how the lift is produced on an 
airfoil-. In order to understand this, the flow around the air- 
foil must be considered. From the Lilienthal effect, we must- 
deduce the conditions of motion of the air, in a way similar to ' 
that described by Knoller and Schmauck, and the meteorological 
cause of such a condition of motion mast then be sought in 
Schnell* s vertical direction of transmission. 

It is now important to get a clear idea of the magnitude of 
the Lilienthal effect on the basis of a possible flow pattern 
and to compare the result with the Knoller-Betz effect alone- 
Corresponding to the method of Prof. W. Hoff (Z.F.M. , Vol. 
13, pp. 276-278), a sinuous wind is taken as the basis, in which, 
however, the ascending flow velocities are greater than the de- 
scending. A streamline must fulfill the condition 

and simultaneously let 

■^o = ^ "^oo n > 1, 
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in wiiicli p'o the greatest angle arid P* the smallest angle 

formed by the streamline with the horizontal at the distance 
from an ordinate 1=0. K = 2 tt/l arid L is the wave length, 
Vq is the maxlraum and Vqo minimum flow velocity. 

Fig. 1 is the flov; diagram which strictly fulfills the con- 
tinuity conditions. It is produced by shifting all the stream- 
line points the same optional small amount K in a direction 
which forms the angle f- with the horizontal. The flow pattern 
can be continued as far as desired perpendicular to the plane of 
the diagram. 



Since 



we have 



^oo ^o 



sin(t - p'o) = 1^ 



or sin(^l- + 3;o) ^ foQ ^ Io_ 
sin(t - P'o) So Vqo 



sin(Tlr + p'o) = J^. 
If we put 

Vq = n Voo* n > 1, 

we then have 

On the horizontal line AB, at the distance I from the ori- 
gin I = 0,,, the flow velocity is 



° sinCt - > 
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and the direction is 

P' = 3»oSin(Xl'), I = r +-AI (3) 

For n = 1 we have the special case of the Knoll er-Betz 
theory. Formulas (l), (3) and (3) then "become 

tan t - T) "t^a-n Po = «» 



(4) 



t = I - 90° 



- = -0^^ - (5) 
= p'o sin(\l) (6) 

since A I = 0, consequently 1=1'. 

For small values of P'q, formula (5) approximates 

and hence the velocity of flow is nearly constant. 

The flow diagram for n = 2 is plotted in Fig. 1. We will 
now endeavor to find what power the same glider can receive from 
the wind in the flow diagram for n = 2 and in the flow diagram 
for n = 1, P o "being 20*^ and the mean wind velocity being the 
same in both cases. For n = 2, we take v - 14 m/s. 

■ 

From formula (1) we obtain 

t = 47.5°, --If - Po = 27.5° 
In Fig. 2 the flow velocities were calculated from formulas 
(2) and (3) multiplied by the cosine of the corresponding angles 
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of flow in texms of I and we obtain the mean horizontal wind 
velocity v^i. 



= i f' V cos p' d I = 9.26 m/s 



o 

We will first consider the case of flight against the wind 
and infinite inertia of the airplane, which is permissible for 
a sufficiently high frequency* Moreover the airfoil constantly 
maintains the angle of attack of least drag. The horizontal ve- 
locity Vp of the airplane will then be such that 

= / d t = 0 
o 

in which T is the duration period, H,^ the drag at the time 
t, also 

and c-^^ is the drag coefficient. 

The present investigation is based on the polar diagram em- 
ployed by Hoff in his article (Z.F»M. 1932, p. 276) and a 16 m 
airfoil* For different values, of v^,, at the time 1;, hence 
at the distance I from the point l = 0 (Fig* l) , the angle of 
attack -p of the wind and its velocity Vj^, before being affect- 
ed by the wing, were determined according to 

+ D _ _vsin 
^^"" P ~ vj, + V cos 0* 

S sin p 
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Puitiiefmorfe^ the c-u values for p at the time t were taken 
from Figi 13^ using the curve o-^ plotted against p. These 
vaiueife weie multiplied by v^ and plotted against t ■ and thus 

was obtained. 

This rendered it possible to plot Hj^^ against Vj, and to 
determine, by the intersection of the connecting curve withtibhe 
abscissa, the value of v for which H__ = 0. In the case un- 
der consideration, it was found to be 

Vj, = 16 m/s 

From Fig. 4 we obtain, for this v^, as the mean value Y^q 
of all the vertical forces, 

^10 = f 0-^^% vjf d t = 506 kg 

i.e., our 16 m® airplane would just be able to maintain horizon- 
tal flight with a full load of 506 kg. The necessary power for 
horizontal flight in still air is approximately ■ ' 



T _ / __F = 506 X 4 / 506 - Q 6? HP 

^ " 75 y p^' 75^ y 16 X 250 " ^'^^ ^' 



in which is the mean flight coefficient, hence the mean of all 
cif/cD^ for the same range of angles of attack, which occur in 
dynamic soaring flight. For our airplane, this flight coeffici- 
ent'.. V^about 250. 

If we now take for our basis a flow pattern in which n = 1 
and the mean wind velocity is the same, i.e., 9.26 m/s, 20^ 
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giving the case of the Betz theory under otherwise the same con- 
ditions. 

We obtain the flow velocities, at the distance I from 
1=0, from 

cos B'o 
V = Vrt — — ^ 
° cos P ' 

and the mean horizontal imgnitude of wind velocity from 

cos B 'n n I o « 

= "^o SSkT' " "^^ COS&0 

The mean horizontal magnitude of wind velocity is therefore 
constant and should be 9.26 ra/s. tt corresponds therefore to the 
magnitude of wind velocity for = 0°* 

The value of Vj,, fot which Hj^^ =0, is found in the same 
way as when n = 2 (Fig. 5). 

We ge-b Vj, = 13 m/s. 

Hio = 0» "^'lo = f-/ Cv^ V ^ * = ^S- 

T _ 368 X 4 / " 568 = r q*^ MP 
^ 75~ y 16 X 250 ^-^^ 

For n = 2, the deducible power is therefore about 62-^ 
greater than for n = 1. This was the most favorable case "under 
the Knoller-Betz theory. In the simplest case the deducible pow- 
er is considerably smaller, but there is a greater percentage in- 
crease for n = 2 as compared with n = 1, due to the camber 
of the airfoil. 

We will now investigate, for the same two flow patterns, the 
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simplest case, in which the airfoil remains in a constant nega- 
tive angle of attack 7. We here assume 7 = - 4°. 

Since I have at my disposal no better profile, which has 
been tested at a sufficiently large negative angle of attache,. I 
will employ the prof ile Rumpler 01 from TB II, p.- 33, though this 
profile is very unfavorable for the case in hand. In Fig. 6, 
c^^ and Cy_^ are plotted against P for this profile* 8 m/s is 
found as the airplane velocity corresponding to the flow pattern 
n = 2, for which the drag vanishes. In Fig. 7, c^ v^^^ and 
^hj^ """r^ Q^® again plotted against t for Vj. = 8 m/s, and we 
have 

^^xo =o/°v, d t = 0 



= ^ / Cv vjf d t = 36 kg. 
J- o 1 



The polar diagram of the profile Ru 01 gives the flight 

coefficient | = cr^/cn^ as 22.5. The mean flight coeffici- 

^ ■'^max 

ent I', therefore, can surely not be over 20, so that, in the 
most unfavorable case, the deducible power will be 



L = /aZF=-36_>L^ / 36_ ^..0.645 HP. 

^ 75 y r 75 y 16 X 20 

in the same manner the result for the flow pattern for 
n. = 1 is found, from Fig. 8 to be 

^10 = 0' ^10 = f o^^ °v, -^R^ ^ ^ = 21.8 kg 
L - \l ^ ^ / !H-fo= 0-203 HP. 
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The increase in the Knoller-Betz effect is therefore, in 
the simplest case, about llbfc for n = 8, as compared with 
n = 1. A considerably greater increase oaJi "be expected vrlth bet- 
ter profiles. 

The relative data obtained Indicate the following: 

1. The best modern gliders (such as the Vampyr and Geheim- 
rath) should, according to their dynamic power requirement of 
about 2.5 HP., be capable of turbulence soaring flight in "the 
most favorable case" of the Knoller-Betz theory. 

Nothing of that kind, however, happened, even in the case 
of airplanes v/hich were supposed to steer automatically according 
to the Betz diagram. From this v/e might conclude either that 
steering according to "the most favorable case" is impossible, or 
that the flow diagram, taken as the basis, corresponds in no way 
to reality. 

2. The "simplest case" affords possible power deductions, 
which, in the evaluation of the flights of the Strolch, Konsul 
and Hawa 6, must have been already very noticeable. If dynamic 
effects were not noticed, however, (i.e., if sinking speed and 
coefficient of glide corresponded, both with wind and in gliding 
flight with no wind, to those theoretically determined for static 
soaring flight, wq might conclude that turbulence is in no way 
capable of improving the flight. In the cases here considered, 
either n or or both values must have been assumed too 
high, or the flow pattern must have been disturbed by other mete- 
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orological causes. 

In opposition to this, however, stands the Lilienthal ef- 
fect, which still declares that the resultant of all the dynamic 
forces, at least in flight against the wind with a velocity 
equal to that of the wind, must be directed upward 3-4?. 

Wc will investigate tho wind itself, with reference to its 
Lilienthal effect, for our case n = 2. We will assume that the 
flow pattern moves with such a velocity against the Lilienthal 
vane that a sufficiently high frequency is produced and that the 
inertia of the experimental arrangement can be regarded as infi- 
nite. The velocity of the flow pattern, in the present investi- 
gation may bo regarded as independent of the wind velocity. It 
simply defines the period or frequency (similar to the propaga- 
tion of water waves ¥/ithout horizontal displacement of the water 
particles) . 

We will resolve the wind of the flow pattern, at the time 
t, into a vertical and a horizontal component. The vertical 
component Vy is 

Vy = V sin P' , ■ 
and the horizontal component 

v^ = V cos P'. 

The continuity is expressed by 

T 

/ V sinp' = 0. 
o 

If we construct the squares of the flow velocities (hence 
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the 16-fold, pressures, in consideration of 

q = P = 0.125) 

and multiply these values "by the sine or cosine of the corre- 
sponding flow inclinations, thus obtaining 

Vy^ sinP' and v^® cos^' 
and plot these values against values of t, we obtain from Fig. 

J / v^ sinp'd t = + 6.0 and i / v^^.^ cosP'd t = 92.0 
" o ^ ^ o 

and thDrefrom the Lilienthal effect through 

tan a = ^ as a = 3° 44' 

Hence, for the flow pattern employed, n and Po^^^e so 
chosen that the magnitude of the Lilienthal effect was limited. 
We must therefore conclude! 

1. That the flow pattern for n = 2, taken as the basis, 
may approximate reality. 

2. That steering according to the Betz diagram for "the 
most favorable case" was not successful. 

3. That the simplest case of the Knoller-Betz theory occurs 
in every soaring flight. The calculated minimum sinking velocity 
for the static soaring flight is represented by the well known 
formula 



D 



max 
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Turbulence of both imperceptible (higher) and perceptible 
frequency is, however, present in every wind. The latter condi- 
tion is demonstrated by the need of steering movements to maintain 
equilibrium. The angle of attack for °L^/°D^ii]a^ gl^^ 
ing flight, be contimiously maintained only in still air. - In a • 
wind, the value of ^L^^jf^^y^ only occuxs incidentally,^ since ' 
then even with a stationary airplane, the angle of attack of the 
wind is subject to continual ciie-nges. The sinking speed must 
therefore be greater in a -v/ind, being approximately "^q ~ J "p^' > 
in which |' is a mean value which denotes, during a period, a max- 
imum value of c-|^®/cjj® attained momentarily at time t. 

According to the quantitative evaluation of the Betz diagram 
in "the most favorable case" by Prof. W. Hoff, we would have, in 
a v/ind, a 25^ greater mininuim sinking speed relative to the wind. 
This did not prove to be the case. Kence: 

4. There is no pure static soaring flight. A portion of the 
flight power is always automatically derived from the turbulence; 
how much is still to be determined. For the present, the djmamic 

power seems to increase the diminution of | = CT^/crf to I* 

^ ■'^max 

and to prevent, by stee'ring motions, the diminution of the coeffi- 
cient of glide. A foothold can be obtained from the data given 
in this article. One can compare the power requirement of the 
best gliders e.g., the Hawa 6, which, with 15 m® area, is 

T - -rr ^ Tns - Q'447 X 145 _ ^ Qf,^ tro 

^stat - 75 ~ 75 0.866 HP . 
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v;ith the results of Figs. 7-S, v;hich give the possible power on 
the basis of winds of 0.6 and 0.3 KP. for 16 m^ area. 

5. There is no use in developing, for the glider profile 

alone, the highest ct Vc-n^ and ct,/ o-n , but only the 

•^max " ^max 

highest mean^values, since soaring flight is attempted only in 
a wind. The vibration amplitudes arc not yet known. Lilienthal 
effect does not produce them since, for the wind here tested, 
n could be greater and S'o smller, without changing the direc- 
tion of the resultant of all the dynamic pressures. The most 
favorable profiles for soaring flight would therefore have to 
bo found by direct measurements of lift and drag on experiment- 
al airfoils in a natural wind. 

Translation by Dwight 11. Miner, 
National Advisory Committee 
for Aeronautics. 



W ^A-. CLJV.^ Tfi-clmloal Memorandum Ho. 324 



Figs.l&2 




H.A.C.A. Technical MemCraiidiM ^6iZ84i 



Figs.3&4 



1.6 
1,4 
1.2 
1.0 
+ .8 
+ i6 
+ .4 
+ .2 
0 









A- 










































































\ 








A 


























































^ 











































































-•4 

-.3 \ 
-.2 

-.1 

-0 

+ .1 
•f..2 



Fig. 3 




'R 



+80 



^ ■f-40 



0 

-40 
j-80 

-120 
-160 



o 



Fig. 4 



N.A.-C.A; Technical Memorandum No. 324 



Figs.5&6 




Fig. 6 R\ampler CI ^=-4° 
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Figs ^7, 8&9- 




Fig. 8 
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